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Abstract 

We present hadron-level predictions from the Monte Carlo generator Cascade and par- 
ton level calculations of open b quark, 6-flavored hadron and inclusive 6-jet production in the 
framework of the fc^-factorization QCD approach for the LHC energies. The unintegrated 
gluon densities in a proton are determined using the CCFM evolution equation and the 
Kimber-Martin-Ryskin (KMR) prescription. Our predictions are compared with the first 
data taken by the CMS and LHCb collaborations at ^/s = 7 TeV. We study the theoretical 
uncertainties of our calculations and investigate the effects coming from parton showers in 
initial and final states. The special role of initial gluon transverse momenta in description 
of the data is pointed out. 

PACS number (s): 12.38.-t, 13.85.-t 



1 Introduction 

Beauty production at high energies is subject of intense studies from both theoretical 
and experimental points of view since events containing b quarks present an important 
background to many of the searches at the LHC. From the theoretical point, the dominant 
production mechanism is believed to be quark pair production through the gluon-gluon 
fusion subprocess and therefore these processes provide an opportunity to test the different 
predictions based on Quantum Chromodynamics (QCD). The present note is motivated 
by the recent measurements [ll|6] of beauty production performed by the CMS and LHCb 
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collaborations at the LHC energy ^/s = 7 TeV. It was observed [2]-S] that the data on 
B^, meson and open h quark production tend to be higher than the MC@NLO [71IB] 
predictions. There are no predictions which describe reasonably well the angular correlations 
between ^-flavored hadrons measured [5] by the CMS collaboration. On the other hand, the 
measurements of transverse momenta and rapidity distributions of beauty hadrons [Ij and 
inclusive 6-jets [6] are reasonably well described by the MC@NLO. 

In the framework of the /c-r-factorization approach of QCD [9], heavy quark production 
has been studied (for previous results see [T014T6] ). In our previous paper ^IQ we have 
obtained a good agreement between the Tevatron data on the open h quarks, bb di-jets, B^ 
and several D mesons (or rather muons from their semileptonic decays) production with the 
predictions coming from fcj^-factorization and we have investigated the role of initial and 
final state parton showers. We have shown that a good description of the specific angular 
correlations between the final-state particles is obtained in Monte Carlo event generator 
Cascade |T7] once the higher order process gg* — ?► gg with subsequent g ^ bb splitting is 
included, which is not discussed here. 

Based on these results, here we give a systematic analysis of the recent CMS and LHCb 
data [IHn] on beauty production in the framework of fc^-factorizatioiJl]. Following to [TB], 
we produce the calculations in two ways: we perform numerical parton-level calculations 
(labeled as LZ) as well as calculations with the full hadron level Monte Carlo event generator 
Cascade and compare both with the measured cross section^. We investigate the influence 
of parton showers in initial and final states for the description of LHC data. Specially we 
concentrate on the angular correlations between the produced 6-fiavored hadrons measured 
by the CMS collaboration [5] which are important for our understanding of production 
dynamics [HHTG]. Finally, we study the different sources of theoretical uncertainties, i.e. 
uncertainties connected with the gluon evolution scheme, heavy quark mass, hard scale of 
partonic subprocess and the heavy quark fragmentation functions. 

The outline of our paper is the following. In Section 2 we recall the basic formulas of the 
/c-r-factorization approach with a brief review of calculation steps. In Section 3 we present 
the numerical results of our calculations and a discussion. Section 4 contains our conclusions. 

2 Theoretical framework 

In the present note we follow the approach described in the earlier publications |14f[TB] . 
For the reader's convenience, we only briefly recall here main points of the theoretical scheme. 
The cross section of heavy quark production in pp collisions at high energies in the kx- 
factorization approach is calculated as a convolution of the off-shell (i.e. /cT-dependent) 
partonic cross section and the unintegrated gluon distributions in a proton. It can be pre- 

iSee also [18] 

^In addition to the comparison of CASCADE predictions with the data in [4] we present here further 
studies. 
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sented in the following form: 



a{pp^QQX) 



I 




(1) 



where ^(x, k|., /x^) is the unintegrated gluon distribution in a proton, \/A{g*g* — )■ is 
the off-shell (i.e. depending on the initial gluon virtualities k^y ^2t) matrix element 
squared and averaged over initial gluon polarizations and colors, and s is the total center-of- 
mass energy. The produced heavy quark Q and anti-quark Q have the transverse momenta 
Pit and p2T and the center-of-mass rapidities yi and ?/2- The initial off-shell gluons have a 
fraction Xi and X2 of the parent protons longitudinal momenta, non-zero transverse momenta 
kiT and (k^T = — fci^ 7^ 0, k^j, = —kl^ 7^ 0) and azimuthal angles 0i and 02. The 
analytic expression for the \M.{g*g* — )■ can be found, for example, in [QlfTBJ. 

The unintegrated gluon distributions in a proton ^(x,ky,/i^) involved in (1) can be 
obtained from the analytical or numerical solutions of the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) [19] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [20j evolution equations. As 
in [12], we have tested a few different sets. First of them, CCFM set AO has been obtained 
in [21] from the CCFM equation where all input parameters have been fitted to describe the 
proton structure function ^2(0:, Q^). Equally good fit was obtained using different values 
for the soft cut and a different value for the width of the intrinsic ky distribution (CCFM 
set BO). Also we will use the unintegrated gluon densities in a proton taken in the Kimber- 
Martin-Ryskin form [22] • The KMR approach is a formalism to construct the unintegrated 
parton distributions from well-known conventional ones. For the input, we have used the 
standard MSTW'2008 (LO) [23] (in LZ calculations) and MRST 99 [21] (in Cascade) sets. 

3 Numerical results 

The unintegrated gluon distributions to be used in the cross section (1) depend on the 
renormalization and factorization scales and ^p. Following to [16], in the numerical 
calculations we set = mq + {v\t + Pir)/^; = ■§ -|- Q|,, where Qr is the transverse 
momentum of the initial off-shell gluon pair, = 1.4 ± 0.1 GeV, mf, = 4.75 ± 0.25 GeV. 
We use the LO formula for the coupling ^^(/il.) with n/ = 4 active quark flavors at Aqcd = 
200 MeV, such that a.(M|) = 0.1232. 

We are in position to present our numerical results. The CMS collaboration has measured 
and meson cross sections in the kinematic range pt{B^) > 5 GeV, \y{B~^)\ < 2.4 [2j 
and pt{B^) > 5 GeV, \y{B^)\ < 2.2 [3j, respectively. The measurements of decay muon 
cross sections have been performed [1] for prifJ^) > 6 GeV and \ri{fi)\ < 2.1. The LHCb col- 
laboration have measured [1] the pseudorapidity distribution of 6-hadrons in forward region 
2 < ri{Hb) < 6, where Ht is any 6-fiavored hadron. In our calculations the fragmentation of 
b quarks into a b hadrons is described with the Peterson fragmentation function ^25j with 
default value of shape parameter = 0.006. To produce muons from b quarks in the LZ cal- 
culations, we first convert b quarks into b hadrons and then simulate their semileptonic decay 
according to the standard electroweak theory. The branching fractions of 6 — B^, b ^ B^, 
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bource 


a{B ) 


a(B ) 




CMS data [fih] 


28.1 ±2.4 ±2.0 ±3.1 


33.2 ±2.5 ±3.5 


1.32 ±0.01 ±0.30 ±0.15 


AO (LZ/Cascade) 


32.7/24.5 


31.4/24.3 


1.31/0.96 


BO (LZ/Cascade) 


23.6/18.8 


22.5/20.4 


0.98/0.72 


KMR (LZ/Cascade) 


16.7/13.1 


15.8/12.4 


0.91/0.59 


MCONLO [2H1] 


19.1 


25.2 


0.95 


Pythia [SHU 


36.2 


49.1 


1.9 



Table 1: The 6- flavored hadron production cross section in pp collisions at ^/s = 7 TeV. 

6 — )■ /i as well as the cascade decay b ^ c ^ fi are taken from [26j . The CMS collaboration 
has presented preliminary data ^ on the inclusive 6-jet production at the i/s = 7 TeV. 
The cross sections have been determined in four 6-jet rapidity regions, namely \y\ < 0.5, 
0.5 < \y\ < 1, 1 < \y\ < 1.5 and 1.5 < \y\ < 2. The 6-jets in the Cascade calculations are 
reconstructed with the anti-fct cone algorithm [27] (using the Fast jet package [251122]) with 
radius R = 0.5. 

The results of our calculations are shown in Figs. 1 — 6 in comparison with the data. 
We obtain a good description of the data when using the CCFM-evolved (namely, AO) gluon 
distribution in LZ calculations. The shape and absolute normalization of measured 6-flavored 
hadron cross sections at forward rapidities are reproduced well (see Fig. 4). The KMR and 
CCFM BO predictions are somewhat below the data. In contrast with b hadron and decay 
muon cross sections, the results for inclusive 6-jet production based on the CCFM and KMR 
gluons are very similar to each other and a reasonable description of the data is obtained by 
all unintegrated gluon distributions under consideration. 

The Cascade predictions tend to lie slightly below the LZ ones and are rather close to 
the MC@NLO calculations [TJI^ (not shown). The observed difference between the LZ and 
Cascade is not surprising and connected with the missing parton shower effects in the LZ 
evaluations. The influence of such effects is clearly demonstrated in Fig. 7, where we show 
separately the results of Cascade calculations without parton shower, with only initial 
state, with only final state and with both initial and final state parton showers. One can 
see that without initial and final state parton showers, the Cascade predictions are very 
close to the LZ ones. The similar situation was pointed out previously [161 Tevatron 
energies. We have checked that the LZ and Cascade predictions coincide at parton level. 

In order to study the dependence of our predictions on the quark-to-hadron fragmentation 
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Source 


a{pp ^b + X^^ + X',p!^>6 GeV, |r/^| < 2.1) 


CCFM set AO 


0.96 /ib 


CCFM set A0+ 


+13% 


CCFM set AO- 


-2% 


mf, = 5.0 GeV 


-7% 


rub = 4.5 GeV 


+6% 




+9% 


Total 


il7% 
=^7% 



Table 2: Systematic uncertainties for beauty total cross section in pp collisions at ^/s = 7 TeV 
obtained with Cascade. 

function, we repeated our calculations with the shifted value of the Peterson shape parameter 
eft = 0.003, which is is often used in the NLO pQCD calculations. Additionally, we have 
applied the non-perturbative fragmentation functions proposed in where the input 

parameters were determined in [2I1E2] by a fit to LEP2 data. The results of our calculations 
are shown in Fig. 8. We find that the predicted cross sections in the considered kinematic 
region are larger for smaller values of the parameter eb or if the fragmentation function 
from [301432] is used. Thus, the LHC data lie within the band of theoretical uncertainties. 

The visible cross sections of 6-fiavored hadrons and 6-decay muons are listed in Table 1 
in comparison with the CMS data pHl]- In Table 2 the systematic uncertainties of our 
calculations are summarized. To estimate the uncertainty coming from the renormalization 
scale fiR, we used the CCFM set A0+ and AO— instead of the default density function AO. 
These two sets represent a variation of the scale used in as in the off-shell matrix element. 
The A0+ stands for a variation of 2 fi^, while set AO— reflects /xr/2. We observe a deviation 
of roughly 13% for set A0+. The uncertainty coming from set AO— is generally smaller and 
negative. The dependence on the 6-quark mass is investigated by variation of 6-quark mass 
of rrib = 4.75 GeV by ±0.25 GeV. The calculated 6-quark cross sections vary by ~ ±6%. 

Now we turn to the investigation of the angular correlations between the produced b 
hadrons. As it was pointed out in [UHlS], such observables are very sensitive to the de- 
tails of the non-collinear gluon evolution. The CMS collaboration ^ has measured the 
6-flavored hadron pair production as a function of the angular separation A0 between the 
two reconstructed b hadrons and variable AR = J (Ar])^ + (A0)2 for three different event 
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scales, characterised by the leading jet transverse momentum pt, namely px > 56 GeV, 
Pt > 84 GeV and pt > 120 GeV. The kinematic range for the measurements is defined by 
the requirements pt{Hi,) > 15 GeV and \ri{Hh)\ < 2 for both of 6-fiavored hadrons. The 
leading jet is required to be within \ri\ < 3.0. Our predictions for and AR distributions 
are shown in Figs. 9 and 10. One can see that none of the calculations fully describes the 
LHC data and therefore there is a room for further studies. Note that the predicted shapes 
of A0 and AR distributions are very different for different unintegrated gluon densities used, 
as it was expected. This is in a contrast to the cross sections as a function of transverse 
momenta or rapidities where all gluon distributions gave a similar behaviour. Note also that 
the measured cross sections at small A0 or AR are significant. Moreover, they exceed the 
ones observed at large angular separation where the two 6-fiavored hadrons are emitted in 
opposite directions. This behavior is reproduced by the KMR gluon density only due to 
different tail at large compared to the CCFM-evolved gluon distributionqfj. The role of 
non-zero gluon transverse momentum k-r is clearly illustrated also in Fig. 11. Here the solid 
histograms correspond to the results obtained according to the master formula (1) and the 
dotted histograms are obtained by using the same formula and without virtualities of the in- 
coming gluons in partonic amplitude. In the last case an additional requirement k^ 2t < A^r 
is applied. One can see that the gluon transverse momentum k-p (both in the hard matrix 
element and in the gluon distribution functions) is important for description of the LHC 
data at low A0 or AR. 

4 Conclusions 

In this note we analyzed the first data on the beauty production in pp collisions at the 
LHC taken by the CMS and LHCb collaborations. Our consideration is based on the /ct- 
factorization approach supplemented with the CCFM-evolved unintegrated gluon densities 
in a proton. The analysis covers the total and differential cross sections of 6-flavored hadrons 
and muons originating from the semileptonic decays of beauty quarks as well as the double 
differential cross sections of inclusive 6-jet production. Using the full hadron-level Monte 
Carlo generator Cascade, we investigated the effects coming from the parton showers in 
initial and final states. Different sources of theoretical uncertainties have been studied. 

Our LZ predictions with the default set of parameters agree with the data on the trans- 
verse momentum and pseudorapidity distributions of 6-flavored hadrons or b quark decay 
muons. The CASCADE predictions tend to slightly underestimate the data at central rapidi- 
ties but the data points still lie within the band of theoretical uncertainties. In this case the 
overall description of the data at a similar level of agreement as in the framework of NLO 
collinear QCD factorization. Special attention was put on the analysis of specific angular 
correlations between the produced 6-fiavored hadrons measured by the CMS collaboration. 
The description of of A(f) and AR distributions in the framework of the /cy-factorization 
require further studies. 

■^A detailed comparison of KMR and CCFM gluon distributions can be found in [33] . 



6 



5 Acknowledgments 



We are grateful to comments from V. Chiochia, G. Dissertori, W. Erdmann, A. Rizzi, 
V. Zhukov and S. Baranov. The authors are very grateful to DESY Directorate for the 
support in the framework of Moscow — DESY project on Monte-Carlo implementation for 
HERA — LHC A.V.L. was supported in part by the grant of president of Russian Federation 
(MK-3977.2011.2). Also this research was supported by the FASl of Russian Federation 
(grant NS-4142.2010.2), FASI state contract 02.740.11.0244, RFBR grant 11-02-01454-a and 
the RMES (grant the Scientific Research on High Energy Physics). 

References 

[1] LHCb CoUaboration, Phys. Lett. B 694, 209 (2010). 

[2] CMS Collaboration, Phys. Rev. Lett. 106, 112001 (2011). 

[3] CMS Collaboration, Phys. Rev. Lett. 106, 252001 (2011). 

[4] CMS Collaboration, JHEP 03, 090 (2011). 

[5] CMS Collaboration, JHEP 1103, 136 (2011). 

[6] V. Chiochia, Nucl. Phys. A 855, 436 (2011). 

[7] S. Frixione and B.R. Webber, JHEP 06, 029 (2002). 

[8] S. Frixione, P. Nason, and B.R. Webber, JHEP 08, 007 (2003). 

[9] L.V. Gribov, E.M. Levin, and M.G. Ryskin, Phys. Rep. 100, 1 (1983); 

E.M. Levin, M.G. Ryskin, Yu.M. Shabelsky, and A.G. Shuvaev, Sov. J. Nucl. Phys. 53, 
657 (1991); 

S. Catani, M. Ciafoloni, and F. Hautmann, Nucl. Phys. B 366, 135 (1991); 
J.C. Collins and R.K. Ellis, Nucl. Phys. B 360, 3 (1991). 

[10] S.P. Baranov and M. Smizanska, Phys. Rev. D 62, 014012 (2000). 

[11] Ph. Hagler, R. Kirschner, A. Schafer, L. Szymanowski, and O.V. Teryaev, Phys. Rev. 
D 62, 071502 (2000). 

[12] M.G. Ryskin, Yu.M. Shabelski, and A.G. Shuvaev, Phys. Atom. Nucl. 64, 1995 (2001); 
Yu.M. Shabelski and A.G. Shuvaev, Phys. Atom. Nucl. 69, 314 (2006). 

[13] N.P. Zotov, A.V. Lipatov, and V.A. Saleev, Phys. Atom. Nucl. 66, 755 (2003). 

[14] S.P. Baranov, N.P. Zotov, and A.V. Lipatov, Phys. Atom. Nucl. 67, 837 (2004). 

[15] A.V. Lipatov, L. L5nnblad, and N.P. Zotov, JHEP 0401, 010 (2004). 

[16] H. Jung, M. Kramer, A.V. Lipatov, and N.P. Zotov, JHEP 1101, 085 (2011). 



7 



[17] H. Jung, Comp. Phys. Comm. 143, 100 (2002); 
H. Jung et al, Eur. Phys. J. C 70, 1237 (2010). 

H. Jung, M. Kramer, A.V. Lipatov, and N.P. Zotov, |arXivillQL4328' [hep-ph]. 

E.A. Kuraev, L.N. Lipatov, and V.S. Fadin, Sov. Phys. JETP 44, 443 (1976); 
E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45, 199 (1977); 
LL Bahtsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 

M. Ciafaloni, Nucl. Phys. B 296, 49 (1988); 

S. Catani, F. Fiorani, and G. Marchesini, Phys. Lett. B 234, 339 (1990); 
S. Catani, F. Fiorani, and G. Marchesini, Nucl Phys. B 336, 18 (1990); 

G. Marchesini, Nucl. Phys. B 445, 49 (1995). 

H. Jung, ii arXivThep-ph/0411287[ 

M.A. Kimber, A.D. Martin and M.G. Ryskin, Phys. Rev. D 63, 114027 (2001); 
G. Watt, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C 31, 73 (2003). 

A.D. Martin, W.J. Stirling, R.S. Thorne, and G. Watt, Eur. Phys. J. C 63, 189 (2009). 

A.D. Martin, R.G. Roberts, W.J. Stirling, and R.S. Thorne, Eur. Phys. J. C 14, 133 
(2000). 

C. Peterson, D. Schlatter, I. Schmitt, and P. Zerwas, Phys. Rev. D 27, 105 (1983). 
K. Hagiwara et al. (PDG Collaboration), Phys. Rev. D 66, 010001 (2002). 
M. Cacciari, G.P. Salam, and G. Soyez, JHEP 0804, 063 (2008). 
M. Cacciari and G.P. Salam, Phys. Lett. B 641, 57 (2006). 



G.S.M. Cacciari and G. Soyez, http://fastjet.fr/ 



M. Cacciari and P. Nason, Phys. Rev. Lett. 89, 122003 (2002); 

M. Cacciari, S. Frixione, M.L. Mangano, P. Nason and G. Ridolfi, JHEP 0407, 033 
(2004). 

M. Cacciari and P. Nason, JHEP 0309, 006 (2003). 

E. Braaten, K.-M. Cheng, S. Fleming and T.C. Yuan, Phys. Rev. D 51, 4819 (1995). 

B. Andersson et al. (Small-x Collaboration), Eur. Phys. J. C 25, 77 (2002); 
J. Andersen et al. (Small-x Collaboration), Eur. Phys. J. C 35, 67 (2004); 
J. Andersen et al. (Small-x Collaboration), Eur. Phys. J. C 48, 53 (2006). 



8 




10 15 20 25 30 

Pt [GeV] 




0) 

O 

F 

m 

t 
a. 
a. 



Q. 

5 10-'^ 
c 
■o 



Cascad e 

• CMS data 

CCFM Set AO 

-• CCFM Set BO 
KMR 



10 15 20 25 30 

[GeV] 



Cascade 



g 15 

m 

T 
a. 

B 10 

>. 

■o 

"B 
■o 



0.5 





• CMS data 




CCFM Set AO ■ 




- - - CCFM Set BO ' 




KMR 






- -t-l 


■ — h—- 





1.5 



B*y 



Figure 1: The transverse momentum and rapidity distributions of meson production at 
the LHC The kinematical cuts apphed are described in the text. The sohd, dashed and 
dotted histograms correspond to the results obtained with the CCFM set AO, BO and KMR 
unintegrated gluon densities. The first column shows the LZ results while the second one 
depicts the Cascade predictions. The experimental data are from CMS [2]. 



9 




^ 10 
0) 

o 

B 

CO 

t 

a 
a 

"B 

■D 



Cascade 



• CMS data 

CCFM Set AO 

-- CCFM Set BO 
KMR 



10 



20 



30 

b" [GeV] 



40 



15 



Q. 

>. 
■a 



10 



Cascade 



0.5 



• CMS data 

CCFM Set AO 

- - CCFM Set BO 
KMR 



1.5 



Figure 2: The transverse momentum and rapidity distributions of meson production at 
the LHC The kinematical cuts apphed are described in the text. The left histograms show 
the LZ numerical results while the right plots depict the Cascade predictions. Notation of 
all histograms is the same as in Fig. 1. The experimental data are from CMS [3]. 
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Figure 3: The transverse momentum and pseudorapidity distributions of muons arising from 
the semileptonic decays of b quarks at the LHC The kinematical cuts apphed are described 
in the text. The left histograms show the LZ numerical results while the right plots depict the 
Cascade predictions. Notation of all histograms is the same as in Fig. 1. The experimental 
data are from CMS [1]. 
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Figure 4: The pseudorapidity distributions of 6-flavored hadrons at the LHC The kinematical 
cuts apphed are described in the text. The left histogram shows the LZ numerical results 
while the right plot depicts the Cascade predictions. Notation of all histograms is the same 
as in Fig. 1. The experimental data are from LHCb pj. 
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Figure 5: The double differential cross sections da/dydpT of inclusive 6-jet production as a 
function of pt in different y regions (LZ predictions). Notation of all histograms is the same 
as in Fig. 1. The experimental data are from CMS [B]. 
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Figure 6: The double differential cross sections da/dydpT of inclusive 6-jet production as a 
function of pt in different y regions (Cascade predictions). Notation of all histograms is 
the same as in Fig. 1. The experimental data are from CMS [B]. 
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Figure 7: Parton shower effects in the pseudorapidity and transverse momentum distributions 
of the 6-quark decay muons. The four hues represent full parton shower (solid line), no parton 
shower (dashed line), initial state parton shower (dotted line) and final state parton shower 
(dashed dotted line). The experimental data are from CMS 
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Figure 8: The dependence of our predictions on the fragmentation scheme (LZ calcula- 
tions). The solid, dashed and dotted histograms correspond to the results obtained using 
the Peterson fragmentation function with = 0.006, = 0.003 and the non-perturbative 
fragmentation functions from [23H2S], respectively. The CCFM-evolved (AO) gluon density 
is applied. The experimental data are from CMS [Mill]- 
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Figure 9: The distributions in azimuthal angle difference between the produced 6-ffavored 
hadrons at the LHC The first column shows the LZ numerical results while the second one 
depicts the Cascade predictions. The kinematical cuts applied are described in the text. 
Notation of all histograms is the same as in Fig. 1. The experimental data are from CMS [5]. 
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Figure 10: The distributions in AR in the 6-flavored hadron production at the LHC. The 
first column shows the LZ numerical results while the second one depicts the Cascade pre- 
dictions. The kinematical cuts applied are described in the text. Notation of all histograms 
is the same as in Fig. 1. The experimental data are from CMS [5]. 
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Figure 11: Importance of non-zero transverse momentum of incoming gluons in open b quark 
production at the LHC The sohd histograms correspond to the results obtained according 
to the master formula (1). The dotted histograms are obtained by using the same formula 
but now we switch off the virtualities of both incoming gluons in partonic amplitude and 
apply an additional requirement k^2T < A^r- have used here the CCFM AO gluon for 
illustration. The experimental data are from CMS [3]. 
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